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Abstract. We present the formalism for computing the critical exponent 
corresponding to the /3-function of the Nambu-Jona-Lasinio model with 
SU{M) X SU{M) continuous chiral symmetry at 0(1/A^^) in a large N 
expansion, where is the number of fermions. We find that the equations 
can only be solved for the case M = 2 and subsequently an analytic expres- 
sion is then derived. This contrasting behaviour between the M = 2 and 
M > 2 cases, which appears first at 0(1/A^), is related to the fact that the 
anomalous dimensions of the bosonic fields are only equivalent for M = 2. 
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1 Introduction. 



Quantum field theories with a four fermi interaction are known to possess 
an interesting range of properties. For example, the two dimensional model, 
introduced by Gross and Neveu in Q , is an asymptotically free theory where 
the fundamental particles of the massless classical theory acquire a mass dy- 
namically in the quantum theory. Therefore, it mimics several of the features 
of more involved four dimensional field theories such as quantum chromo- 
dynamics and so has been used as a simple laboratory for testing out ideas 
which are harder to examine there. More recently, four dimensional mod- 
els possessing a four fermi interaction have received intense study following 
the observation of Nambu, Q, that such models could provide a realistic 
alternative to the Higgs boson for generating a mass for the particles seen 
in nature, |^ 

One technique which can be used to analyse models with a four fermi 
interaction is the large expansion, Q, where N is the number of fermions. 
For instance, computing the effective potential of the two dimensional model 
in the saddle point approximation when is large one observes that the 
perturbative vacuum of the theory is not the correct one but that where 
the vacuum expectation value of a bosonic auxiliary field is non-zero which 
in turn generates the fermion mass. This behaviour is also preserved in 
the three dimensional model and it is important to note this property is 
non-perturbative and therefore would never be accessed in a perturbative 
approach. However, whilst the large N expansion has been powerful in 
revealing such structure the conventional method possesses one particular 
drawback in that it is virtually impossible to calculate beyond the leading 
order. This is primarily due to the nature of the (dynamically generated) 
propagators of the bosonic ipip bound states which exist in the model and 
are non- fundamental in structure, Thus computing Feynman integrals 
with such massive propagators becomes impossible beyond leading order. 

Clearly there is a need to probe such models to higher orders in and 
to have a technique which interpolates between the spacetime dimensions 
already discussed. One such method currently does exist which achieves 
these aims is the critical point self consistency method of Vasil'ev et al 
which was developed for the 0{N) bosonic a model in |5|, The beauty 
of that approach is that one calculates within the theory at the fixed point 
which is defined as the non-trivial zero of the d-dimensional /^-function. 
At this point the theory is conformal and therefore massless which allows 
one to overcome the difficulty of computing previously intractable integrals 
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beyond the leading order. Moreover, one solves the model by determining 
the critical exponents of the fields and Green's functions |p. These can then 
be related to the perturbative renormalization group functions through an 
examination of the renormalization group equation in the critical region. 
Examples of this approach can be found in Further, the /5- function of any 
theory carries all the important information on the evolution of the coupling 
constant with the renormalization scale, and it too can be computed order 
by order in \/N by considering the appropriate critical exponent, Q. As 
these are determined as functions of the spacetime dimension one can deduce 
information on models simultaneously in several different dimensions. 

In this paper, we complete the examination of the four fermi models 
at 0(l/iV^) by computing the /3- function exponent of the Nambu-Jona- 
Lasinio model with a non-abelian continuous global chiral symmetry. Pre- 
viously various exponents had been calculated in the 0{N) Gross Neveu 
model in [^, ^ 0] using the original large approach. The break- 
through toO(l/iV2) was achieved in with the computation of the fermion 
anomalous dimension 77 and latterly the mass, [13, 14 1, vertex, |ll3| , 
and /3-function exponents, [0, ^ have been determined as well as rj at 
0(1/A^^), p^, p^ . Several of these exponents were subsequently determined 
to 0{1/N'^), Ijl^, for the generalization of the 0{N) model, which possesses 
a discrete chiral symmetry, to the case where it has aC/(l)xC/(l)or SU{M) 
X SU{M) continuous global chiral symmetry, [Q, Several leading order 
results had earlier been presented in |11, ll^ using the saddle point ap- 



proach. By computing the exponent A, where 2A = — P'{gc), at 0{1/N'^) 
in the SU{M) x SU{M) case here the model can then be said to have 
been solved thermodynamically since knowledge of two independent expo- 
nents means that the remaining ones can be deduced through hyperscaling 
laws which have been checked at leading order in [11| and their consistency 
merely reflects the renormalizability of the model. Several additional mo- 
tivations for considering the SU{M) x SU (M) model in particular include 
the provision of analytic results for the three dimensional model which will 
provide key estimates for numerical simulations of the model on the lat- 
tice. Recently, estimates for several exponents have been determined using 
Monte Carlo methods in the 0{N) Gross Neveu model for relatively small 
values of N in [20| and they are in agreement with 0{1/N'^) results, |15]. 
Further the provision of expressions valid in d-dimensions as a function of 
is important for demonstrating the equivalence of various models. For 
example, it is known that the (4 — e)-dimensional Yukawa model and the 
(2 + e)-dimensional Gross Neveu model lie in the same universality classes 
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which has been estabUshed by examining the e-expansions of the perturba- 
tive renormahzation group functions to as many orders as they are presently 
known. In our case the equivalence of the (2 + e)-dimensional non-abelian 
Nambu-Jona-Lasinio or chiral Gross Neveu model will probably be with the 
generalized (4 — e)-dimensional Gell-Mann-Levy a model introduced in [^] 
and currently of interest due to its relation to hadronic physics. (The origi- 
nal model possessed an SU (2) x SU{2) chiral symmetry and by generalized 
we mean its natural extension to the case of SU{M) x SU{M).) The pro- 
vision of critical exponents to 0{1/N'^) which we carry out here will be the 
first step in such a proof. 

Finally, we mention that there is some uncertainty in the literature con- 
cerning the behaviour of the SU{M) x SU{M) model for M > 2 and for 



M = 2, 1 11]. It has been suggested in |ll|] that the chiral symmetry is only 
realised in the case M = 2 and not for higher M, IQ. So far in the large 
self consistency approach there has been no indication of distinct behaviour 
between either case. This is normally observed in the breakdown of the 
computation of the /3-function exponent at some order in large N which has 
been discussed in the context of other models in |22|. At leading order we 
were able in [17| to compute r/, A and the vertex anomalous dimensions, Xi 
at 0{1/N) and r] and x ^-t 0(l/iV^) without complications for all M. So it 
will be interesting to see if it is possible to deduce A at 0{1/N'^) for all M. 

The paper is organised as follows. In section 2 we introduce the details 
of the model we will examine at criticality and illustrate the method by 
computing r/ at 0{1/N). The leading order analysis is continued in section 
3 where the calculation of A at 0{1/N) is detailed in preparation for the 
0{1/N'^) analysis which is presented in section 4. We conclude the paper 
by solving the master equation in section 5 where we discuss the results for 
M = 2 and M > 2 separately. 



2 Preliminaries. 

To begin with we describe the model we are interested and introduce the 
notation and formalism we will use in solving it. Its lagrangian can be 
written in several ways but we choose to use the auxiliary field version 
since it involves three point vertices which are essential to the technique of 



uniqueness used to compute the Feynman graphs, [^. Thus we take 



L = iV5^^^^^^ + aV;^V^ + i7rV^A?j7V'^-^(a2 + ^»2) (2.1) 

25 
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where ■0*^ is the fermion field with l<i<N,l<I, J< M, 1 < a < 
(M^ — 1) and g is the coupHng constant. Eliminating the auxiliary bosonic 
fields a and tt"" from (2.1) would yield the four fermi interaction explicitly 
and according to (2.1) possesses an SU{M) x SU{M) continuous global 
chiral symmetry. The matrices A"j are the generalized Pauli matrices of the 
group SU (M), where [A*^, A^] = 2i/"'"^A'^ and /"'"^ are its structure constants, 
and we normalize them in as general a way as possible by taking 

Tr(A"A'') = 4T{R)6''^ (2.2) 

The Casimirs which will appear throughout the paper are 

A-^A" = AC2{R)I , f'""^f'"^ = 02(0)6"'' (2.3) 

where for SU{M), T{R) = \, 02{R) = (M^ - 1)/2M and 02{0) = M. 
Also, in order to carry out intermediate checks on our calculation we will 
retain 02{R) and 02(0) in their general form to allow us to take the limits 
to either the case with U{1) x U{1) chiral symmetry, (C2(i?) = T{R) = 
J, 02(0) = 0, M = 1) or the 0{N) model with discrete chiral symmetry, 
{02{R) = 02(0) = 0) and compare with results we already know. 

We now introduce the basics to enable us to solve (2.1) to 0{1/N'^) in 
the critical point self consistency approach. First, the fundamental idea is 
to examine the scaling behaviour of the fields of (2.1) in the critical region 
defined as the non-trivial zero of the /3-function in d-dimensions, |^, ^ 
Since the fields are massless at this fixed point where there is a conformal 
symmetry the structure of the Green's functions take a particularly simple 
form. For instance, in coordinate space they are 

^(^) ~ ' ^(^) ~ ' ^(^) ~ (2.4) 

as — >■ where we use the same letter to denote the propagator and A, B 
and O are x-independent amplitudes. By analogy with ideas in statistical 
mechanics the properties of the model in the critical neighbourhood are de- 
scribed totally by the critical exponents q, /3 and 7 of each of the individual 
fields, 0. By the universality principle they are functions only of the space- 
time dimension, d = 2jjL, and any internal parameters of the theory. For 
our case, these will be and M through the appearance of the Casimirs of 
(2.3). Moreover, the exponents can be related to the appropriate renormal- 
ization group functions which are ordinarily calculated perturbatively as a 
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series in g, by examining the renormalization group equation at criticality, 
[0]. For instance, if we examine the canonical dimension of each field in the 
action which is a dimensionless quantity then we can define their anomalous 
dimension as follows, [p^] , 

a = fi + ^rj , 13 = l-r]-Xa , j = I - rj - Xw (2.5) 

where rj is the fermion anomalous dimension which is related through the 
critical renormalization group equation to the fermion wave function renor- 
malization constant. The other exponents, Xa and Xtd are the anomalous 
dimensions of the respective 3- vertices in (2.1) and these can also be related 
to the analogous renormalization constants. 

In the critical point approach we discuss throughout this paper, the 
exponents are expanded order by order in where is large, and the 

exponent rji, say, is calculated at the ith order where r] = J2'^i Vi/^^ and r/j 
= r]i{fi,M). The method to do this has been discussed extensively before, 
[p2i, but for completeness sake and for preparation to calculating the (3- 



function exponent A we review the calculation of ryi, [17|. This is achieved 
by examining the skeleton Dyson equations with dressed propagators in the 
critical region, truncated to the appropriate order in which for the 

moment will be leading order. The relevant Feynman graphs for this are 
illustrated in fig. 1, where the quantities cr~^ and tt~^ are the 2-point 
functions of the fields. Their asymptotic scaling forms in the critical region 
in coordinate space can be easily deduced from the forms of the propagators 
(2.4) by mapping (2.4) first to momentum space performing the inversion 
there before inverting the map to coordinate space, ||5[ This is facilitated 
by use of the Fourier transform 

1 a(a) f e^'^^ ^^^^ 



(a;2)" 22"7r'^ (A;2)M-a 
where, for simplicity, a(a) = r(/i — a)/T{a). The functions are then 
N r{a - 1):^ _^ p{[3) 

^ A(x2)2m-«+1 ' ~ S(x2)2m-/3 



C(x2)2m-7 



clS X where 
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To represent the Dyson equations of fig. 1 in the critical region, one merely 
substitutes (2.4) in the lines of the one loop graphs to obtain 

= r{a - 1) + z + 4C2iR)y (2-9) 
= p{(3) + 2zNM (2.10) 
= p{j) + 8T{R)Ny (2.11) 

where z = A^B and y = A^C and we have used the relations (2.2) and (2.3). 
In (2.9)-(2.11) one has three unknowns, z, y and iji through the appearance 
of T{fj, — a) in (2.9). Thus eliminating zi and yi after expanding to leading 
order in 1/A^ one obtains 

c2{Ry 



^1 = TT 



2 



1 



where fji 



T{R) 

i)/[r(/i + i)r(/.)r(i-/i)r(/. 
fj,r^{fi)m Mzi 



47r2MM 



yi 



AT{R) 



(2.12) 
1)] and consequently 
(2.13) 



which will be required later. Whilst (2.12) is a result for all M and agrees 
in the limits mentioned earlier with js], ^ |ll| and |]l^], it also checks with 
a leading order calculation in [11 1 for M = 2, using canonical techniques. 
However, as has been mentioned this latter method is not powerful enough 
to probe to 0{l/N'^). 

The vertex anomalous dimensions of (2.5) have also been computed to 
the same order in 1/A^ by examining the scaling behaviour of the 3-point 
functions in an analogous way at criticality by extending the earlier work of 
[p4| , 25] and we record the results in preparation for computing A, as 



Xal 



XttI 



2(^ - 1) 
2{fx - 1) 



1 

M 



C2iR) 
T{R) 

C2{R) 1 
T{R) M 



C2{G) 



2T{R) 



(2.14) 
(2.15) 



Again each agrees in the appropriate limit with earlier results and for M 
2, Xo-i agrees with pT], 



3 /3-function exponent. 

Having indicated the simplicity of the method to deduce the exponents of all 
the fields at leading order, we now present the extensions to the formalism 
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to deduce Ai which will serve as the foundation for calculating A2 where A 
= /U — 1 + J2iZi Ai/A^*. To achieve this one considers the corrections to the 
asymptotic scaling forms (2.4) by including O(x^) terms, |5|, ^, |l2|, 



V'(^) ~ + (3-1) 

B 



^(^) ~ + (3.2) 

<x) ~ ^[l + C'{x')'] (3.3) 

Additional higher order corrections, which will involve other exponents such 
as the specific heat, could also be computed by including the appropriate 
forms. However, as there exist hyperscaling laws relating such exponents, 
by doing this one would merely succeed in demonstrating their consistency 
which is not in doubt, [||, The quantities A', B' and C are the respec- 
tive amplitudes associated with each correction and the analogous 2-point 
functions can be deduced in a similar way to (2.7) to obtain, [p^ ], 

V'-'w - :4^§)v^[i-^'-(«-i)(^')'] (3-4) 

^-H^) - B^jy,-, [^-B'm{^'f] (3.5) 

^"'(^) - cX^^[l-C^'«(7)(^')'] (3.6) 

as ^ where 

a(a-/i + A)a(a-A) q{a)a{a - ^) 

= 7 r-r-^ > s{a) = — — — (3.7) 

a{a — fj,)a[a) [a — n + \)[a — \) 

With these corrections one can now reconsider the Dyson equations in 
the critical region. However, when one is dealing with a theory where the 
fundamental field is fermionic in this large N approach, as we have here, 
several higher order graphs have to be included in the a and vr Dyson equa- 
tions, |22|, which are illustrated in fig. 2. The reason for this is quite simple 
and can be seen, for instance, by examining the A'^-dependence of each term 
in the cj-consistency equation. First, we have from figs 1 and 2, in the critical 
region 

= B'q{p){x'^)^]+2zM[l + 2A'{x'^)^] 
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- z^M[U + {x^)\A'U,a + B'Uib)] 

+ 4yzMC2{R)[U + {x^)^{A'UiA + C'Uic)] (3.8) 

where the subscripts A, B and C denote the respective insertion of (x^)"^ on 
a ■0) cr or TT hne. As in ||6|, ^ the equation decouples into two pieces one of 
which is relevant for the computation of 172, which we ignore here, and the 
other which involves the (x^)'*' terms, since each is of differing dimension in 



2 • 

X , le 



= AzA' - B' 



+ 4C2iR)yzUw (3.9) 



NM 

In (3.9), we have neglected the graphs where there is an insertion on the 
fermion line of the 2-loop graphs of fig. 2 as these do not contribute at 
0{1/N). Further, the expressions Hia, and Hie denote the value of 
the integral which is 0(1), without symmetry factors. Now if one substitutes 
the leading order values for a, (3 and 7 into (3.9) to examine the location 
of in each term, it is easy to see that both terms of the coefficient of 
B' are of the same order. Therefore, neglecting the graph Hib would omit 
a contribution and lead to an erroneous result and so it must be included. 
Similarly, the vr consistency equation gives 



= lQT{R)yA' + 4T{R)yzkiBB' 



C 



P^MjI + 8T{R)y\2C2{R) - C2{G))Aic 



(3.10) 



Finally, in the ip equation one needs only to consider the one loop graphs 
of fig. 1 since there is no reordering as there is in (3.9) and (3.10) which can 
be seen by inspection. Thus 

= A'[-r{a-l)s{a-l)+z + 4yC2{R)]+zB' + AyC2{R)C' (3.11) 

To proceed one forms a 3 x 3 matrix with A', B' and C as the basis vectors 
and sets its determinant to zero for consistency whence an expression for Ai 
will emerge since it is the only unknown. For the moment the matrix is 

/ r{a - l)s{a - 1) z C2{R)y \ 

4z PMm + ,2ji> -C2{R)yzn' 

A,, pMqM I y^n'{2C2(R)-C2{G)) 

\ y^^^ iqt(r)n 2 / 

(3.12) 

where the explicit calculation of the 2-loop graphs gives Hib = Hic = Ai^ 
= Aic = n' = 27r2^/[r2(^)(^-l)2], m, 11. Again it is easy to see by exam- 
ining the A-dependence in each element of (3.12) the necessity of including 
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the graphs of fig. 2. To solve for Ai by setting the determinant of (3.12) to 
zero one makes several row and column transformations. Subtracting row 
two from row three and then column three from column two one obtains the 
following matrix whose determinant is set to zero 



r(a — l)s{a — 1) 

-4 2qi(5)-U'[ 



^1^-^^ MC2{R) \ 



MC2{R)z 
2T(R) 
MC2{R)zU' 
2T{R) 
Mzn'{4C2iR)-C2iG)) 

(3.13) 

where we have used (2.13). Thus setting (3.12) to zero yields Ai as 



V 







T{R) J 
, _ MC2{R) \ 

^ r(i?) ) 




K7) 



Ai 



{2^Ji - 1)7/1 



1 c^jR) 

M T{R) 



(3.14) 



This agrees with previous results of the 0{N) model [O] and the U{1) 



U{1) case |11, 26] in the appropriate limits as well as the SU{2) x SU{2) 
calculation of However, it is worth noting that an alternative second 
solution might appear to emerge from the lower right element of (3.13). We 
discard this as it is not consistent with previous results and, moreover, the 
row and column transformations we have made do not induce any contribu- 
tion from the ip equations which appear as the upper row in (3.12). Indeed 
this type of factorization of the final determinant has been observed in other 
contexts, [22|, where one had also to ignore a spurious solution. Also, we 



remark that together with (2.12), (2.14) and (2.15) we now have a complete 
set of exponents for the SU{M) x SU{M) model at leading order. 

We conclude this section by stating that it is now possible to proceed 
beyond (3.14) and attempt to calculate the 0{1/N'^) corrections, which has 
been achieved recently for the 0{N) model ||l^, ^ and U{1) x U{1) case 
^ and use can be made of results calculated there. Something which is 



necessary for this is the 0{1/N ) corrections to r], Xa and x-k since these 
will appear in the 1/A^ expansion of the functions of (3.4)-(3.6). They have 
been deduced in ||l^ by considering the scaling behaviour of the higher order 
graphs in the Dyson equations of the 2 and 3 point functions with dressed 



propagators and we record the results here. We have, [17| 



V2 = -r 



+ 



1 I C2iR)y 

M T{R) J 



1 

1^2il 



(3.15) 




T{R)J 



C2{G)C2{R) 
2T'^{R) 



^(/^) + 



2(^ - 1) 
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(the sign of the coefficient of the C2 (G) term was incorrectly given in [T^ ) 



Xa2 



and 



X-K2 



+ 



+ 



+ 



ml 



+ 



4(/.-l)2 
liC2{R)C2{G) 



+ 



+ 




(3.16) 



4(^ - 1)2 
X (2^- 



^(/^) + 



1 



C2{R) 



fC2{R) 



T{R) J V T{R) 
I^C2{G) fC2{R) 2_ _ C2{G) \ 



1 G2{G) \ 
M 2T{R) ) 



2T{R) \T{R) 
3^l^i (C2{R) 



M 2T{R) J 

1 G2(G) Y 
M 2T{R) ) 



G2{G) 



4(/.-l)3 VT(i?) 

L (Chm _ }_ 

16(/x-1)3m[ \T{R) M) T{R) 

( 1 G2{G) f C2{R) 

' I M2 2T{R)M \T{R) 



2ixm 



+ — 



1 

16 



1)^ 



C2{G) y 

2T{R)) 



'(C2{R) 1 Y 



G2{G) _ 3G2{G)G2{R) C|(G) \ 

Mr(i?) 2r2(i?) 2r2(i2) i 



16 Vm r(ii) 

G2(G) / C2(^) 
AT{R) V r(i?) 

+ ^r3e(,) 



g2(G) Y 
2T{R)) 



(2/^2 -5^ + 4) ^ 3 ^ C2(i^) 



+ 



3 

M 

(2^^ 

(m - 1)^ 



( C2{R) 
\T{R) 

3)- 




VM2 T{R)M 
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1 C2{G) \ 
M 2T{R) ) 



C2{R) 
T{R) 



C2{G) \ 
4T{R) J 



(3.17) 



where ^'(/i) = ij{2fi - 1) - ipil) + ip{2 - fi) - and G(/i) = ip'{fi) - 

tp'{l) where 'ip{x) is the logarithmic derivative of the F-function. Further, 
from the r/2 consistency equations 



Z2 



87r2MM 



X + 



1 



2/2 



/"r2(/x)7?2 



1 



C2{R) 
T{R) 



+ 



327r2Mr(ii) 



1 



C2iR) 



1 C2{R) 

{fi-l)KM T{R) 

( 1 C2{R)\ 
1)2 \M T{R) ) 



T{R) 



+ 



M 



(3.18) 



( C2{R) C2{G) \ 

ifi-l) V T{R) M 2T{R) J 
f G2{R) 1 C2{G) \ 
{fi-l)J (/X - 1)2 r(i?) M 2T{R)) 



(3.19) 



4 Beyond 0(l/iV). 



The main effort required to go beyond (3.14) and determine A2 involves 
computing the higher order Feynman graphs which appear in the Dyson 
equations of each field. First, though we need to derive the formal cor- 
rections to each critical representation of the Dyson equations which is a 
non-trivial exercise due to the appearance of divergent graphs which have 
a vertex subgraph and finite three and four loop graphs. The renormaliza- 
tion techniques required to handle the infinities have been discussed in other 
places, |l2|, but for completeness we detail the procedure here and con- 
centrate on the "0 equation since there are fewer corrections to be considered 
compared to the a and vr fields because of the extra graphs which arise. The 
higher order graphs with dressed propagators for the ip Dyson equations are 
illustrated in fig. 3. Including them in (2.9) and (3.11) one obtains 







r(a - 1)[1 - A's{a - l)ix'^)^] + z[l + {A' + B'){x^)^]u^{x^)'>^'^+'^ 
+ 4C2{R)y[l + {A' + C")(x2)^]t;2(x2)x-+^ 

+ z2[S + (A'Sa + B'Sb)(x2)^](x2)2^-+2^ (4.1) 

- 8yzC2iR)[^ + (A'Sa + B'J^B + C"Sc)(x2)^](x2)>^'^+>^-+2^ 

+ 16y^C2iR)[C2iR) - iC2(G)][S + (A'^a + C'Sc)(x2)^](x2)2x.+2A 

Unlike at leading order the powers of do not cancel here and, moreover, 
the graphs S, S^, and Sc are divergent with each having simple poles 



12 



in A. This infinitesimal quantity is a regularization introduced to control 
such infinities by shifting the exponents of a and 7rby/3— >/? — A,7^7 — 
A, [|6|, |l2|. If we formally define the A-finite part of the two loop corrections 
to be etc the simple poles in the (x^)^ sector of (4.1) are absorbed by 
choosing the respective vertex counterterms u and v in such a way that 
(4.1) is A-finite at 0{1/N'^). After this the Inx^ terms which remain and 
which prevent one approaching the — > region are removed by choosing 
Xai and XttI appropriately. It turns out that the choices already found in 
the analysis of the 3-point functions, (2.14) and (2.15), achieve this which 
is a consistency check on our renormalization. Decoupling the equation as 
before yields 

= A'[-r{a-l)s{a-l) + z + iyC2iR) 

+ ^'a{z^ - 8yzC2{R) + l6yHC2{R) - iC2(G)))] 

+ B'[z + zi:'B{z-8yC2iR))] 

+ C'[4C2{R)y + 8yC2{R)T.'c{y{2C2{R) - C2{G)) - z)] (4.2) 

which is the 0{1/N'^) correction to (3.11). However, examining the A^- 
dependence of the terms in the coefficient of A' , is 0(1/A^) relative 
to r(a — l)s{a — 1) and can therefore be neglected in the overall 3x3 
determinant. Moreover, the explicit calculation of and reveals that 
both are zero ^% |l^] and therefore (4.2) simplifies to 

= A'[-r{a-l)s{a-l) + z + AyC2{R)] + zB' + AC2{R)yC' (4.3) 

where of course now the 0(1/A^^) terms of z and y are needed. We have 
detailed the renormalization of the ip equation since it has a simpler structure 
than that of the a and vr equations which are complicated by the higher 
order graphs which have to be included. The relevant graphs, with dressed 
propagators, for the a equation are illustrated in fig. 4 where the label 
beside each graph will be used in the following. We have only displayed 
the distinct topologies in fig. 4 and the solid internal lines, without dots, 
denote either a a or tt field. As at leading order we need only consider [x^)'^ 
insertions on the bosonic lines, since within the final 3x3 matrix graphs 
with insertions on the fermionic lines will contribute to A3 only. The dot is 
intended to indicate the location of the {x^)'^ insertion when one substitutes 
the asymptotic scaling forms of the propagators (3.1)-(3.3) in the graphs. 
The correction graphs for the vr equations are formally equivalent to those 
of fig. 4. It is clear from the structure of 11251 and 1135 that they are A- 
divergent due to the presence of a vertex subgraph. Within the consistency 
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equation the same renormalization procedure and vertex counterterms we 
discussed for the ip equation removes the simple poles in A and it is therefore 
a straightforward matter to write down the finite 0{1/N'^) correction to (3.9) 
in the critical region as 



N 



= zMA'[4- zUia]- B' 
+ C'[AC2iR)Myznic-Ii-C2] 
where IIb2 and Ilc2 are formally defined to be 



+ Z^MUiB + Ub2 



(4.4) 



ns2 
nc2 



2 (nasi + U2B2) + + U4B 
2{2U5Bi + II5B2) - meBi + 2UeB2] 

= 2(n2ci + n2C2) + nsc + 

- 2(2n5Cl + n5C2) - 4(n6C71 + 2n6C2) 



(4.5) 
(4.6) 



with the A-finite part only understood as contributing. The factors of z 
and y as well as group factors are contained within each formal term of 
(4.5) and (4.6) and will be computed explicitly later. Whilst Hi^ and Hie 
contributed at leading order it is important to remember that each has an 
0{1/N) correction which will appear in A2. Further, explicit calculation has 



shown that n'. 



lA 



0, 



The formal corrections to the tt equation are also 



now straightforward to write down and we find 







l6T{R)yA' + B'[AT{R)yzAiB - Ab2] 



P{l)(l{l) 
N 



+ Sy''T{R){2C2{R) - C2{G))K^c + Ac2 



(4.7) 



where the higher order corrections, Kb2 and Ac2, have the same formal 
definitions as (4.5) and (4.6). Thus, with (4.2), (4.4) and (4.7) one can now 
form the corrected version of the 3x3 matrix (3.12) and deduce the master 
equation which yields A2. All that remains is the explicit evaluation of the 
higher order graphs. 

These fall into two classes. First, we need to consider the 0{\/N) cor- 
rections to HiB, Hic, KiB and Aic since, for example, the fermionic lines 
of the 2-loop graphs each have exponents ^ + \ri which therefore means 
they have to be expanded in powers of 1/A^. An algorithm to do this was 



presented in |27| and the analogous graph has also been computed in the 
0{N) model, [|l2[ Whilst the graphs are of a similar structure the ex- 
ponents of the internal bosonic lines differ in each case and we therefore had 
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to recompute Hib = ^ic and Ai^ = Aic explicitly. Using the techniques 
of 1 27] we obtained 



n 



IB 



+ 



(^-i)2r2(^) 

3//(/i - l)r?i 
2N 



1 



l)N \M 
1 



l_ ^ C2{R) 



1 



1 



T{R) 
C2{R) 



1 , C2{R) 
M T{R) 2{fi -1) \M T{R) 



A 



2^T^^' 



IB 



+ 



(^-l)2r2(^) 



1 



2N 



J_ , C2{R) 
M 



l)N \M 
1 



l_ ^ C2{R) 



T{R) 



1 ( C2{R) C2{G) \ 
T{R) 2{i2 - 1) V T{R) M 2T{R) ) 



(4.9) 



The second class of graphs are those of fig. 4. Since the integral structure 
of each is completely the same as the graphs of the 0{N) model we needed 
only to compute the SU{M) factors multiplying each which arise from the 
appearance of A"7^ at various vertices of the graphs. Making use of the 
relations, which have been discussed in |2q|, 



4T{R) 



abc \ c 



abc \ c 



M 



which defines the totally symmetric tensor d"'^'^ in SU{M) and 



^apq^bpq 

fjO-PQ jbrp jcrq 
^apq^brp^crq 



4C2iR) - - C2iG) 



M 



!:ab 



C2(G') ^afec 



4^2 (i?) 



8T{R) 3C2(G) 



M 



'abc 



(4.10) 

(4.11) 
(4.12) 
(4.13) 



we managed to compute 1152) ^C2i ^B2 and Kc2 relative to the 0{N) model 
graphs, which are denoted by a superscript ° in the following. Making use 
of the relation between yi and zi we have 



Ub2 = 2z'M-'{^ + ^^ 
\M T{R) 



[11252 + 
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- ziM(2n?Bi + n?B2 + 2ngBi) 



nc2 = 



+ 2/m2 



Ji ~ T^) + ~ 4ziMngB2] (4.14) 



T{R) 



1 

M 



C2(i?) , C2(G)\ 



r(it;) 



+ 



n9 



2B1 



2T{R) ) 

M T{R) 



n§B (4.15) 



+ 



1_ CM _ C2{G)\ o 
M T{R) 2T{R)) 



Ab2 



+ 



■ 2M2/ 



1 

M 



^^^^^^)n§. (4.16) 



r(i?) J \M T{R) 2T{R) 



Ac2 = 2z^M 



\T{R) 2T{R))\T{R) M 2T{R) J 



+ 



+ 



+ 



l\T{R) 2T{R))\T{R) M 'H\j 

( C2{R) (Chim _ C2{G)\ 2\ o 
\MT{R) \T{R) 2T{R)J J 

( C2{R) _ C2{G)_\ /^CR) _ J_ _ C2{G)_\ o 
\T{R) 2T{R)J\T{R) M 2T{R)) 
( C2{R) _ C^G) \ (l_ qm _ C2{G)\ o 
\T{R) 2T{R) J \M ^ T{R) 2T{R)) 



C2{G) \ /J_ 

2T{R)) \M T{R) 



- ziM 



+ 2 



+ 



[R) 2T(R)J\M T[R) 2T(R) J 

// J_ C2{R) ( C2{R) _ J_ _ C2{G)y Cl{G) \ 

y^W^ MT{R) \T{R) M 4T{R)J 16T^{R)J 

X (2n§si + u2b2) 

2(1. f J_ + ^kCR) _ ^M") 



+ 4 



M \M T{R) 2T{R)) 

( C2{R) C2{G) \ f C2{R) G2{G) \\ 

V T{R) M AT{R) J V T{R) M 2T{R) J ) 

(C2{R) 1 C2{G)\(C2{R) G2{G)\ ^ \1 

\T{R) M 2T{R)J\T{R) AT{R) J ^^vJ ^' ^ 



16 



The explicit values of the basic graphs of the 0{N) model are 

2 



^2B1 
11252 

^5B2 



^6B2 



+ 



+ 



+ 



+ 







- i)3r4(^)A 








- i)^r4(^) 








- i)3r4(;u)A 







1 - A(/i- 1) 36 + 



36 + 



A(^-l) 



1 



36 + 



(2/i - 3)7r6A'o(2/i - 2) 
2(^-l)5(;u-2)r3(;.) 



(4.18) 
(4.19) 
(4.20) 
(4.21) 



1 



+ 



66 - $ - ^2 + 



(2//-3)(/i-l) (2/x-3)(/i-l) 
7r6Ma(2/z-2) [(2/^-3) ^ 



+ 



+ 



(/x-2) 
1 



(^-i)5r3(^) 

(3// -4)^ 
(/.-l)(M-2)2 ' (/.-2)2 
7r6^a(2;u - 2) 

7r<^^a2(2/i - 2) 



2(/i-l)2 (2/X-3) 
2(^-2)^ , (/i-2) 

1 



+ 



+ 



1 



2(/x-l)2 
27r6A'a2(2/i - 2) 
(/^-l)'^(^-2)2 
5 (2/i-l)^' 



(/x-1) 2(/i-l)2 



(2/x - 3)7r6Ma(2^ - 2) 
(/i-l)5(/.-2) 
^' 1 



(4.22) 



(4.23) 



+ 



£ ^ 3(/i-l) 
2" T ~ (2/i - 3) 
(/. - 2)^^ 



6 + 



2(^-1) 4(/i-l)2 (2^-3)(^-l) 



(/i-2) 



+ 



2(2/i-3)(M-l)2 (2/i-3)(/.-l) 



■K^^'a^{2ix - 2) 
(/z - 1)7(^ - 2) 



(4.24) 



where $(^) = ip'{2fj, — 1) — ■0'(2 — /x) — '(/''(/Li) + •0'(1)- It is worth recalling 
that in the U{1) x U{1) case, [^], only three of the above were relevant for 
A2 due to the high degree of symmetry present in that model. 
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5 Discussion. 



We are now in a position to substitute for all the quantities in the overall 3 
X 3 matrix and set its determinant to zero. If one follows the same formal 
transformations which were made at leading order to simplify the evaluation 
we notice that the factorization of the matrix into two solutions, one of which 
is irrelevant, does not occur at 0{1/N'^). For instance, the (13) element is, 
up to factors, 

/x2r2(/i)r/i \2C2{R) 2 C2(G)- 



4T{R)y - Mz 



8{fi - l)7r2MiV2 



T{R) M 2T{R) 



(.-1) <^-" 

which is zero at leading order but not at 0{1/N'^). Similarly the (23) element 
also vanishes at 0{1/N) but not at the subsequent order for all M. This is 
a rather unfortunate situation especially given the fact that it was possible 
to write down 0(1/A^^) expressions for -q, Xa and Xir and A at 0{1/N) for 
all M for this model. The absence of a solution can, however, be related 
to what we will term as the non-realisation of the chiral symmetry of the 
model. For instance, it is easy to see that when M = 2 in the complicated 
expressions (3.16) and (3.17), Xa- = X-n- to 0{l/N'^), but this simple relation 
does not hold for M > 2, ||. In the U{1) x ^7(1) model, 0, X<7 is also 
equivalent to x-k at 0{1/N'^) which indicated that the chiral symmetry is 
preserved in that instance. Now if one examines (5.1) it is zero at M = 2 and 
no other (positive) value of M for all fi since the group factor is proportional 
to (M^ -4). Whilst the different behaviour of the M = 2 and M > 2 models 
has been suggested in [11] where the non-realisation of the chiral symmetry 
was postulated this is as far as we are aware the first observation of it in an 
explicit exponent calculation. 

We now restrict attention to M = 2 as this is the only case we can 
obtain A2 in all dimensions. (Of course, it would be possible to derive an 
expression for A2 for M > 2 as the solution of a quadratic equation, but 
the complicated and non-linear nature of such an exponent would not be in 
keeping with the structure of exponents found in other models. Indeed a 
similar situation has also arisen before in |22] when the non-factorization of 
the consistency determinant indicated that that exponent was not a sensible 
quantity to compute.) It is worth recording the explicit values of the various 
exponents when M = 2 in order to solve (4.2), (4.4) and (4.7) as, fl% , 

m = m (5-2) 
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Xu2 



m 

X-K2 



XttI 



2(/i - 1) + 2/i + 



3/x 



1) 4(m-1)2 



2(/i-l) 
3^(/i- l)e + 2(/i-2)^' 



(5.3) 
(5.4) 



Moreover, 

nc2 

Ai32 
AC2 



8(/x-l)2 
(5;u- l)(2/i^ -5/i + 4) 

{2^x - 1)771 



8z3[2(n?B2 + - 2zi(2n?Bi + n?B2 + 2n?Bi) 

- (n§Bi + n§B - 8zin?B2)] 

24^3 [n^Bi + n§B - 2zi(2n?Bi + n?B2 - 2yUbi)\ 
8z3[n§Bi + - 2zi(2n^Bi + n?B2 - 2ngBi)] 
8z3[n§Bi + 2n§B2 + n§B + 2n2B 

- 8zi(2n?Bi + n?B2-ngB2)] 



(5.5) 
(5.6) 



from which it is easy to verify that 

n_B2 + nc2 — A_B2 



Ac2 







(5.7) 
(5.8) 
(5.9) 

(5.10) 



(5.11) 



which imphes that the (23) element of the transformed determinant is zero 
at 0(l/iV^) as at leading order and as in the U{\) x U{\) case, |2^, which 
again reinforces our point of view on the importance of the chiral symmetry 
at 0(1/A^^). Hence the formal correction to (3.13) is 



[r(a- l)s(a- 1) -4z] 



N 



32^2 (5.12) 



With the values given earlier, some tedious algebra leads to 

^(6/^2 _ 12^ + 5) ;u(2/x-3)($ + ^'2) 



A2 



2(/i - 2)2r?i 



+ 



+ 



3/i(2/i - 3)(2^ - l)G(^) 5 



1 



8(/i - 2) 
37/i 



+ 



1 



2/i 2;u2 8(/i - 2)2 



4(/. - 1) 



23 

4 ^ 16(// - 



+ 4/x2 + - 



55 



8 16(/x - 2) 
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,(f<-2)2 (p-2) 
Also, from (5.13) one can restrict to three dimensions to find 



(5.13) 



1 16 , 8(277r^ + 125) 
^ - 2 " 3^ + 277r4iV2 ^^'^^^ 

Indeed, it is worth comparing (5.14) to the 0(1/A^^) versions of A in other 
four fermi models to observe a reassuring similarity in their structure. For 
instance, in the 0{N) model, |14, |l5|] , 



1 16 32(27vr^ + 632) 
" 2 3^+ 27;^^iV^ 



and in the U{1) x C/(l) case, |26|] , 



1 16 16(277r2 + 472) 
^ - 2 " 3^ + 27^4iV2 (5-16) 

Thus, numerically the exponents differ first at 0{1/N'^) in each model. 

To conclude with we remark that we have discovered that the first oc- 
curence of distinct behaviour between the M = 2 and M > 2 models arises 
at 0{1/N'^) which could therefore not be noticed in a canonical large N ap- 
proach. This feature is intimately related to the equivalence or otherwise of 
the anomalous dimensions of the a and tt fields, and therefore the realisation 
of a chiral symmetry. However, the deeper implications that this important 
observation has in relation to four dimensional field theories is beyond the 
scope of the present paper but does deserve consideration. 
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Figure Captions. 

Fig. 1. Leading order skeleton Dyson equations with dressed propagators. 

Fig. 2. Additional graphs for Ai. 

Fig. 3. Higher order corrections for ip equation. 

Fig. 4. 0{1/N'^) contributions to a equation. 
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